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We develop a renormalization method for calculating the electronic structure of single and double 
quantum dots under intense ac fields. The nanostructures are emulated by lattice models with a 
clear continuum limit of the effective-mass and single-particle approximations. The coupling to the 
ac field is treated non-perturbatively by means of the Floquet Hamiltonian. The renormalization 
approach allows the study of dressed states of the nanoscopic system with realistic geometries as 
well arbitrary strong ac fields. We give examples of a single quantum dot, emphasizing the analysis 
of the effective-mass limit for lattice models, and double-dot structures, where we discuss the limit 
of the well used two-level approximation. 



I. INTRODUCTION 



The renormalization method, as a tool for studying the 
electronic structure of solids, has attracted attention al- 
ready twenty years ago with, the application to disordered 
one-dimensional systemsEn It is a powerful approach for 
systems whose complexity hinders the direct diagonaliza- 
tion of the corresponding Hamiltonian. At that time its 
complete usefulness could be hardly realized due to the 
difficulty in applying the method simultaneously in more 
than one spatial direction. More recently its suitability 
has been revealed for studying strongly anisotropic solids 
that started to deserve growing interest, lika.semiconduc- 
tor superlattices and conducting polymers. El 

In the present work the renormalization method is ap- 
plied to a different problem of increasing interest with 
the same formal structure: the dressed electronic spectra 
of quantum dot systems under intense ac fields. 

In semiconductor technology, the continuous improve- 
ment on growth techniques opened possibilities to de- 
sign systems that are in the quantum limit in all spa- 
tial directions, establishing an important and pioneering 
branch of what is nowadays called nanoscience. One of 
these devices, a double-quantum dot is a candidate for 
the elementary unit of the foreseen quantum computing. 
Two coupled .quantum dots afterwards called as artifi- 
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Mesoscopic systems, like quantum dots and quantum 
dot arrays, are described by lattice models, treated in a 
tight-binding framework, emulating the, continuum limit 
of the effective-mass approximation. 13 Lattice models 
have been used mainly in the context of disorder effects 
on electronic and transport properties of two-dimensional 
systems and quantum billiards.E_l On the other hand, 
there are reports in the literature on using lattice models 
in another situation, where they constitute also an ap- 
propriate approach, namely the simulation of arrays of 
quantum dots and antidots in the presence of a magnetic 



fieldHI Realistic description of quantum dot nanostruc- 
tures in the independent particle approximation inter- 
acting with intense ac fields becomes possible within the 
proposed renormalization-decimation procedure.cl 

The coupling to an ac field is non-perturbatively in- 
cluded using the Floquet method, by means of a pro- 
cedure introduced by Shirley.Ea The Floquet method has 
been widely used for the non-perturbative study of the in- 
teraction of atomic, molecular and semiconductor nanos- 
tructures systems with a strong ac electric field. The 
time-independent infinite matrix Hamiltonian obtained 
after the application of the Floquet-Fourier transforma- 
tion over the time-dependent Schrodinger equation, de- 
scribes entirely these processes without any further ad- 
hoc assumptions. 

Therefore, the effect of an intense ac field on the elec- 
tronic spectra of a nanostructure, like a quantum dot, is 
well described by an infinite Floquet matrix. However, 
due to the large vector basis used to define the system 
in a lattice-model tight-binding approximation, the nec- 
essary convergence criteria make the eigenvalue calcula- 
tions practically impossible beyond the perturbative field 
intensity range. Moreover, we observe that the infinite 
matrix Hamiltonian structure resembles the structure of 
a linear chain matrix, where the "energy sites" are en- 
ergy sub-matrices which corresponds to the Hamiltonian 
of the system plus or minus an associated multiple of 
the photon energy and the "hopping parameter" is the 
ac field coupling the sub-matrices. These relationships 
allow us to develop an interesting and promising renor- 
malization approach whereby the actual dimension of the 
system to be calculated is reduced to the one of the lat- 
tice model for the bare system. Besides that, since the 
quantity calculated is the density of states, one gets a 
step further than by the direct diagonalization of the Flo- 
quet Hamiltonian, which provides only the quasi-energy 
spectra without the spectral modulation as a function of 
the field strength. The hierarchy of these quasi-energy 
spectra related to different photon replicas will also be 
discussed. 

The paper is presented as follows. A brief introduction 
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to the Floquet Hamiltonian is given in Sec. II, followed 
by a comprehensive description of the renormalization 
method applied to the system. In Sec. Ill, numerical cal- 
culations for single and double quantum dots are shown 
and discussed. We first focus on the validity of a lattice 
model for a quantum dot in the presence of intense ac 
fields. The second part is centered on the ac field de- 
pendence of the bonding and antibonding states of an 
artificial molecule, i.e., a double-quantum dot structure. 
In Sec. IV we present our conclusions. 



II. THEORY 
A. The Floquet Hamiltonian 

The lower part of the energy spectrum of a mesoscopic 
system, like a quantum dot or quantum-dot array de- 
scribed in the framework of the effective-mass approxima- 
tion, will be evaluated here using a tight-binding model 
for a square lattice of s-like orbitals, considering only 
nearest-neighbors interaction and a hopping parameter 
defined by V = -h 2 /{2 m*a , ), where m* is the effective 
mass and a is the host lattice parameter. 

The applied ac fields are parallel to one of the square 
sides. Hence, our model is described by the Hamiltonian 
H = H + H int where 
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with (Ti lt i 3 = \h,l2 > and aj j =< h,h\- The atomic 
energies will be taken constant, ei lt i 2 = 4\V\, for all sites. 
F and w are the ac field amplitude and frequency respec- 
tively, and e is the electron charge. The treatment of 
the time-dependent problem is based on Floquet states 
\l±, I2, m > where l\, I2 are the site indexes and m is the 
photori-j|ndex. We follow the procedure developed by 
ShirleyEa which consists in a Fourier-Floquet transfor- 
mation of the time-dependent Hamiltonian into a time- 
independent infinite matrix which must be truncated. 
The matrix elements are 
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where Fi = ^eaF. The dimension of the matrix is 
L\ x L 2 (2M + 1), where L x and L 2 are the maximum 
number of atomic sites, while M is the maximum photon 
index. We choose M in order to satisfy a convergence 
condition: symmetric spectra relative to the edges of the 
quasi Brillouin zones (QBZs). The first QBZ is spanned 
in the range —hui/2 < £ < Tiu/2 . 

The truncated Floquet matrix resembles a linear chain 
matrix in the form 
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where E m = £ — (S + mhio) is a L\ x L2 matrix, with £ 
being the quasi-energy spectrum and S, is the L\ x L 2 
system matrix. For other hand, the L\ x L 2 diagonal 
matrix J- represents the coupling of the system with the 
ac electric field and is given by 



T = FiliS Vih di 2 i 2 



(5) 



A time-independent Schriidinger equation defined in 
terms of Green functions asE3 



{£ - H)G = 1 



(6) 



can be associated to the Floquet matrix, as we will show 
next. 



B. Renormalization Method 

Equation (||) can be projected to the photon Hilbert 
spaces < n\ and \m >, resulting in 



< n\(£ - H)\k > G km = 5 n 
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Expanding this equation, considering £ — > £ + in, where 
77 — ► 0, eliminating the intermediate index sites, making 
some substitutions, and repeating the decimation process 
successively we arrive to the general expression, for any 
£ > 1 order decimation of the Eq. ([?]): 
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The Green function Goo after the first decimation is 
jiven by 



Goo — 
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where 
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The density of states of Floquet is expressed as 

p(£ + i V ) = lim ( - - Im [ Tr Goo ] ) • (10) 



III. ARTIFICIAL ATOMS AND MOLECULES 

A. Quantum dot under intense ac field 

As a starting point we calculate the quasi-density of 
states of quantum dots as a function of the ac field in- 
tensity. We choose a square and a rounded geometry 
to represent the quantum dot as shown in Fig. |l| The 
square quantum dot consists of an array of 6 x 6 atomic 
sites, Fig. |l|(a). A more realistic geometry is given in 
Fig.[I](b), which will be used for the study of double-dot 
systems in this work. It should be noticed that the area 
of this geometry (37 sites) is comparable with the square 



one if we use the same tight-binding parameters given be- 
low. Such lattice models exhibit a particle-hole symme- 
try in the electronic structure and are usually thought as 
simple, although useful, approximations for superlattices 
or arrays of quantum dots, where each quantum well or 
quantum dot is represented by a lattice site, respectively. 
This extreme lattice limit has been used for studying 
qualitativeUt.the effect of intense ac fields on superlattice 
minibands.c3 On the other hand, lattice models may be 
useful in calculating the lower part of electronic systems 
well described by the effective mass approximation. In 
the present work, the tight-binding hopping parameter is 
chosen in order to emulate the electron effective mass for 
the GaAs bottom of the conduction band, to* = 0.067m . 



Since V 



-0.142 eV for a model 
lattice parameter of a — 20 A. This leads to a quan- 
tum dot with a lateral width of L = 120 A, an order of 
magnitude lower than typical dimensions of actual quan- 
tum dots fabricated by lithographic methods. We intend 
to illustrate the method with the present calculations. 
Nevertheless, all conclusions can be traced back to larger 
systems by an appropriate scaling of the energy parame- 
ters. 

The first point to be addressed here is establishing the 
limit between the extreme-lattice and effective-mass lim- 
its. This is achieved by following the evolution of the 
quasi-density of states as a function of the ac field inten- 
sity with frequencies of the order of the electronic band 
width. In FigJ| we show the contour plot of the density 
of states as a function of the field intensity for the bot- 
tom half of the first QBZ. We recall that, for the chosen 
hopping parameter, the band width for a square array 
of single atomic s-like showing particle-hole symmetry 
is given by AE = 8\V\ = 1.136 eV. Here the field fre- 
quency is Tilj = 1.0 eV. The square quantum dot, even 
though a textbook example, is very useful to understand 
the breakdown of the effective-mass limit. At zero field 
intensity the energy spectrum of a 2D square potential 
well with L = 120 A - a quantum dot - is clearly identi- 
fied, as well as the breakdown of degeneracies due to the 
applied ac field. The ac Stark shifts lead to crossings at 
eaF/Tuj w 2.4, corresponding to a dynamic localization 
along the chains perpendicular to the field direction. This 
is the limit where the host lattice effects are already pre- 
dominant, i.e., the discrete basis of atomic orbitals does 
not emulate the effective-mass approximation anymore 
and the length scale is given by the host lattice param- 
eter a and not by the lateral width L of the quantum 
dot. Indeed, the model simulates in this limit a square 
array of quantum dots, each one represented by a sin- 
gle site, resembling the spectra that would be expected 
for coupled chains, where each chain mimics the dynamic 
localization in superlattice minibands. 

This result for a square array is a good guide for charac- 
terizing the density of states as a function of the field in- 
tensity of a quantum dot with lower symmetry, Fig 
as can be seen in the equivalent plot shown in Fig 



(b), 
|(b), 

also for Tiui = 1.0 eV. The dynamic localization at the ex- 
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treme lattice limit is less denned due to the fact that the 
atomic site chains perpendicular to the field direction are 
not equivalent as in the square array. Furthermore some 
of the degeneracies are already broken in the absence of 
the ac field. 

One clear advantage of the method can already be seen 
in Fig. ||. By diagonalizing the Floquet Hamiltonian, the 
quasi-energy spectra are depicted in the so called QBZs, 
each of them reproducing the spectrum with the energy 
shifted by integer multiples of the photon energy. The 
overlap of these photon replicas makes the interpreta- 
tion of the spectrum rather cumbersome, specially for 
strong overlaps, which are unavoidable for low frequen- 
cies. The present renormalization approach indicates 
that the quasi-energy spectrum, modulated by the field 
dependent density of states for higher or lower photon 
replicas becomes relevant only with increasing field in- 
tensity. This effect is verified at the bottom of Figs. ||(a) 
and @(b). The spectrum replica lowered by one photon 
energy shows a negligible contribution at low field. On 
the other hand, in both figures we see no significant mod- 
ulation of the density of states in the low field limit within 
the entire band of the depicted zero photon replica, due 
to the high frequency considered. 

The main interest, however, is the effective-mass limit, 
i.e., the energy bottom of each photon replica at low 
fields in the scale of Fig. ||, as well as low frequencies, 
that would couple only these few low-energy states. The 
scaling of these quantities shows the suitability of the 
present method, as will be seen in the following discus- 
sion on double quantum dots. 



B. Double Quantum dot 

Our double-quantum dot system is based on the quan- 
tum dot shown in Fig. 0(b). The coupling between dots 
is of free choice and we consider a simple connection of 
both dots by the same hopping parameters considered 



so far, as illustrated in Fig. 



Such configuration is an 



example of strong inter-dot coupling. The corresponding 
quasi-density of states plot as a function of the field in- 
tensity for the same high field frequency, Tvjj — 1.0 eV, 
as in Fig. || is shown in Fig. ||(a) . The quasi-density of 
states spectrum is very similar with the appearance of 
the expected energy splitting due to the coupling among 
quantum dot states. The structure of the energy split- 
ting may be rather complex, considering the coupling 
of initially degenerate states in each quantum dot. A 
covalent-like binding, with a splitting between a bonding 
and antibonding states, is well defined for the inter-dot 
coupling between the lowest state in each quantum dot. 
Therefore, from now on we will focus exclusively on the 
energy bottom of the spectra, the continuum low-field in- 
tensity limit, in order to analyze the lowest pair of split 
double-quantum dot states. 

In Fig. ||(b) we have the quasi-density of states for the 
lowest pair of molecular states. Both bonding and an- 



tibonding states shift rigidly upwards in energy. Here a 
modulation of the density of states is seen. By increas- 
ing the field intensity, the density of states diminishes, 
with increasing contribution of higher and lower photon 
replica (not shown). Since the frequency of the field is 
very high hu> = 1 eV, all states of the system are mixed 
by the field and no typical two-level behavior is observed. 
The rigid energy shift may be seen already as a lattice 
effect, since for high field intensities, the dynamic local- 
ization for the host lattice is observed, Fig.[|(a). 

A clear covalent picture is revealed for much lower 
frequencies, in the range of Tilo w 10 meV, which is of 
the order of the tunnel splitting between the lowest pair 
of states in the bare "molecule" for the chosen param- 
eters: A sp iit = 7.1 meV. For this frequency range the 
coupling to higher molecular states is negligible, since 
the third molecular state is about 50 meV above the an- 
tibonding state. Lattice effects are also absent since the 
dynamic localization effects on the host lattice are rele- 
vant for field frequencies of the order of the entire spec- 
tral width, which is two orders of magnitude larger than 
the energy scale of interest given by the tunnel splitting. 
The present calculations are exact, since we are in an 
independent particle approximation, for an artificial 
molecule. Within this parameter range, the lowest pair 
of states of the double-quantum dot behave as a two-level 
system, as will be discussed next. We are going to switch 
the field frequency from above to below the tunnel split- 
ting energy. The quasi-density of states plots will also 
be compared to the corresponding quasi-energy spectra 
of the equivalent two level system, obtained by diagonal- 
izing the Floquet Hamiltonian. 

In Fig. 5 the evolution of the artificial bonding 
and antibonding states as a function of the field inten- 
sity is shown for a frequency fuv = 10 meV > A sp ii t . In 
Fig. 5(a) we see the quasi-density of states around the 
main "photon replica" , which shows the highest intensity 
for low fields. The next important branches of the quasi- 
density of states are the bonding state plus one photon 
and the antibonding minus one photon. With increasing 
the intensity, the splitting between the zero-photon bond- 
ing and antibonding states diminishes down to a crossing 
at eaF/hui « 0.3. This resembles the dynamic localiza- 
tion regime for superlattice minibands.ES The branches 
of the quasi-density of states in Fig. 5(a) can be mapped 
on the quasi-energy spectrum of a two-level system, em- 
ulated by two atomic sites with fitted tight-binding pa- 
rameters. This spectrum is a function of a field inten- 
sity defined by edF/fru;, where d is the distance between 
the two effective atomic sites. The correspondence be- 
tween Figs. 5(a) and 5(b) is satisfied by properly scaling 
d w 7a. Having in mind the double quantum dot of Fig. 
3, do — 7a is the distance between the centers of the dots, 
i.e. "the bond length of the molecule" . 

The equivalent situation, for a frequency given by 
htu = 7 meV, therefore near the resonance situation 
A S piit ~ hoj, is shown in Fig. 6. Here we clearly see 
that the bonding and antibonding states evolve in Rabi 
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sidebands!^ This behavior can also be mapped on an 
effective two-level system, Fig. 6(b). The effective dis- 
tance between the two effective atomic sites is close to 
the "molecular bond length" do = 7a, like for frequencies 
higher than the tunnel splitting, Fig. 5. 

The quasi-resonant case may be analyzed directly from 
the Rabi frequency at the crossing of the sidebands in 
Fig. 6(a). At the crossing Hu>r = A sp ut, where ujr = 
dnF/Ti with dn being the dipole matrix element. Since 
Aspiit = fko, one has that eaF w O.ihojR, resulting in a 
dipole matrix element d# rs 3.33a and 2dn = 6.6a, in 
good agreement with the "molecular length" do = 7a. 

In Fig. 7 we show the situation for a field frequency 
Tuo = 5 meV lower than the bare tunnelling splitting. 
The zero-photon molecular states show an ac Stark shift, 
Fig. 7(a). The dressed tunnel splitting increases with 
field intensity up to the first important anticrossing at 
eaF/hu = 1. Now the mapping on an effective two-level 
system, Fig. 7(b), occurs for an effective distance d « 
3.5, nearly half the nominal molecular bond length. This 
shrinking of the effective distance is, however, compatible 
with the increasing of the tunnel splitting. 

IV. CONCLUSIONS 

The present results demonstrate the usefulness of the 
renormalization-decimation method applied to the Flo- 



quet Hamiltonian. Such a procedure makes possible 
the calculation of quasi-density of states of realistic ge- 
ometries for nanostructures under intense ac fields de- 
scribed by lattice models. The matrix dimension in the 
renormalization method is L\ x L-2, while the direct di- 
agonalization has to handle with matrix dimension of 
(2M + 1) x (Lx x L 2 )- This is of paramount importance 
if we must consider M ~ 100 for L\ j— > 20 in order 
to simulate realistic mesoscopic systems £3 The examples 
shown here for double quantum dots, although heuris- 
tic, point out interesting dependence of the dressed elec- 
tronic structure on the field frequency. A natural exten- 
sion of the work is the analysis of the local quasi-density 
of states, as well as the study of asymmetric double- 
quantum dots, a situation for which the density of states 
reveals actual tunnelling properties, in order to discuss 
transport measurements in such systems. 
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FIG. 1: Quantum dot geometries: (a) square and (b) rounded 
square dots, where a is the host lattice parameter. 



FIG. 2: Density of states for the (a) square and (b) rounded 
dots as a function of the field intensity. The field frequency is 
Tiu) = 1.0 eV. Only the bottom half of the first QZB is shown 
(see text). 



FIG. 3: Double quantum-dot geometry based on rounded 
square dots, where do is the "bond length" . 



FIG. 4: Density of states of the double quantum-dot shown in 
Fig.H as a function of field intensity at high frequency Tu) = 1 
eV, (a) the bottom half of the first QBZ and (b) the bonding 
and antibonding states. 



FIG. 5: (a) Density of states of a double quantum-dot as a 
function of the field intensity for a frequency, hui — 10 meV; 
(b) spectrum of a equivalent two-level system, simulated with 
a dipole distance d w do, in the dynamic localization regime. 



FIG. 6: (a) Same as Fig. ^| for hcj — 7 meV, and (b) same as 
Fig. n near the Rabi resonance. 



FIG. 7: (a) same as Fig. |5por huj = 5 meV, and (b) same as 
Fig. B in the ac Stark regime. 
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